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Abstract 

Real emission corrections to Higgs production via gluon fusion, at order a^, lead to a 
Higgs plus two-jet final state. We present the calculation of these scattering amplitudes, 
as induced by top-quark triangle-, box- and pentagon-loop diagrams. These diagrams 
are evaluated analytically for arbitrary top mass mt. We study the renormalization 
and factorization scale-dependence of the resulting H + 2 jet cross section, and discuss 
phenomenologically important distributions at the LHC. The gluon fusion results are 
compared to expectations for weak-boson fusion cross sections. 



1 Introduction 

Gluon fusion and weak-boson fusion are expected to be the most copious sources of Higgs 
bosons in pp-coUisions at the Large Hadron Colhder (LHC) at CERN. Beyond representing 
the most promising discovery processes [jl|, H, these two production modes are also expected 
to provide a weahh of information on Higgs couphngs to gauge bosons and fermions [§] . The 
extraction of Higgs boson couphngs, in particular, requires precise predictions of production 
cross sections. 

Next-to- leading order (NLO) QCD corrections to the inclusive gluon-fusion cross section 
are known to be large, leading to a if-factor close to two 0. Because the lowest order 
process is loop induced, a full NNLO calculation would entail a three-loop evaluation, which 
presently is not feasible. In the intermediate Higgs mass range, which is favored by elec- 
troweak precision data 0], the Higgs boson mass uih is small compared to the top-quark pair 
threshold and the large limit promises to be an adequate approximation. Consequently, 
present efforts on a NNLO calculation of the inclusive gluon-fusion cross section concentrate 
on the rrit -^oo limit, in which the task reduces to an effective two-loop calculation ^j. In 
order to assess the validity of this approximation, gluon-fusion cross-section calculations, 
which include all finite rrit corrections, are needed. Of particular interest are phase space 
regions where one or several of the kinematical invariants are of the order of, or exceed, the 
top-quark mass, i.e. regions of large Higgs boson or jet transverse momenta, or regions where 
dijet invariant masses become large. For larger Higgs boson masses, top-mass corrections 
become important and a full calculation of if + 2 jet production is needed. 

A key component of the program to measure Higgs boson couplings at the LHC is the 
weak-boson fusion (WBF) process, qq^qqH via t-channel W or Z exchange, characterized by 
two forward quark jets 0]. QCD radiative corrections to WBF are known to be small [0 and, 
hence, this process promises small systematic errors, if -|- 2 jet production via gluon fusion, 
while part of the inclusive Higgs signal, constitutes a background when trying to isolate the 
HWW and HZZ couplings responsible for the WBF process. A precise description of this 
background is needed in order to separate the two major sources of if + 2 jet events: one 
needs to find characteristic distributions which distinguish the weak boson fusion process 
from gluon fusion. One such feature is the typical large invariant mass of the two quark 
jets in WBF. A priori, this large kinematic invariant, m'jj ^ 4m^, invalidates the heavy top 
approximation and requires a full evaluation of all top-mass effects. We will find, however, 
that even in this phase-space region the large rrit limit works extremely well, provided that 
jet transverse momenta remain small compared to rrit. 

In a previous letter ^ we presented first results of our evaluation of the real-emission cor- 
rections to gluon fusion which lead to H + 2 parton final states, at order a^. The contributing 
subprocesses include quark-quark scattering which involves top-quark triangles, quark-gluon 
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scattering processes which are mediated by top-quark triangles and boxes, and gluon scatter- 
ing which requires pentagon diagrams in addition. The purpose of this paper is to provide 
details of our calculation and to give a more complete discussion of its phenomenological 
implications. In Section ^, we start with a brief overview of the calculation. Full expressions 
for the quark-quark and the quark-gluon scattering amplitudes are given in Section ^ Ex- 
pressions for the gg^ggH amplitudes, which were obtained by symbolic manipulation, are 
too long to be given exphcitly. Instead we describe the details of the calculational procedure 
in Section p.3| . The matrix elements for all subprocesses have been checked both analytically 
and numerically. The most important of these tests are described in Section ^. We then turn 
to numerical results, in particular to implications for LHC phenomenology. In Section ^, we 
first compare overall H + 2 iet cross sections from weak-boson fusion and from gluon fusion 
and determine the subprocess decomposition of the latter. QCD uncertainties are assessed 
via a discussion of the scale dependence (renormalization and factorization) of our results. 
We investigate various distributions, searching for characteristic differences between gluon 
fusion and WBF. Our final conclusions are given in Section p. 

A number of technical details are collected in the Appendixes. Scalar integrals, in par- 
ticular the evaluation of scalar five-point functions, are discussed in Appendix ^ Ap- 
pendix ^ gives useful relations among Passarino-Veltman Cij and Dij functions. Finally, in 
Appendixes 0, 0, and |^, we provide expressions for the color decomposition and the tensor 
integrals encountered in triangle, box and pentagon graphs. 

2 Outline of the calculation 

The production of a Higgs boson in association with two jets, at order a^, can proceed 
via the subprocesses 

qq-yqqH, qQ^qQH, qg^qgH, gg^ggH, (2.1) 

and all crossing-related processes. Here the first two entries denote scattering of identical 
and non-identical quark flavors. In Fig. |I| we have collected a few representative Feynman 
diagrams which contribute to subprocesses with four, two and zero external quarks. In our 
calculation, the top quark is treated as massive, but we neglect all other quark masses, so 
that the Higgs boson only couples via a top-quark loop. Typically we have a ggH coupling 
through a triangle loop (Fig. (a)), a gggH coupling mediated by a box loop (Fig. |I] (b)) 
and a ggggH coupling which is induced through a pentagon loop (Fig. Q (c)). The number 
and type of Feynman diagrams can be easily built from the simpler dijet QCD production 
processes at leading order. One needs to insert the Higgs-gluon "vertices" into the tree-level 
diagrams for 2 — 2 QCD parton scattering in all possible ways. 

In the following counting, we exploit Furry's theorem, i.e. we are counting as one the two 
charge-conjugation related diagrams where the loop momentum is running clockwise and 
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Figure 1: Examples of Feynman graphs contributing to H + 2 jet production via 
gluon fusion. 

counter-clockwise. This halves the number of diagrams. In addition, the crossed processes 
are not listed as extra diagrams, but are included in the final results. Three distinct classes 
of processes need to be considered. 

1. qq — > qqH and qQ — > qQH There are only 2 diagrams obtained from the insertion 
of a triangle loop into the tree-level diagrams for qq qq. One of them is depicted in 
Fig. (a), while the other is obtained by interchanging the two identical final quarks. 
In the case of qQ qQH., where Q is a different flavor, there is only one diagram, i.e. 
Fig. g (a). 

2. qg — >qgH At tree level, there are 3 diagrams contributing to the process qg^qg'. 
one with a three-gluon vertex and two Compton-like ones. Inserting a triangle loop 
into every gluon line, we have a total of 7 different diagrams. In addition, we can insert 
a box loop into the diagram with the three-gluon vertex, in 3 different ways: the 3! 
permutations of the 3 gluons are reduced to 3 graphs by using Furry's theorem. In 
total we have 10 different diagrams for the qg qgH scattering amplitude. 

3- qg — ^ggH Four diagrams contribute to the tree-level scattering process gg^gg'. a 
four-gluon vertex diagram and 3 diagrams with two three-gluon vertices each. Inserting 
a triangle loop in any of the gluonic legs gives rise to 19 different diagrams. The 
insertion of the box loop in the 3 diagrams with three-gluon vertices yields another 18 
diagrams. Finally, there are 12 pentagon diagrams (corresponding to 4! permutations 
of the external gluons, divided by 2, according to Furry's theorem). 

The amplitudes for these processes are ultraviolet and infrared flnite in D = A dimensions. 
Nevertheless, we kept D arbitrary in several parts of our computation because some functions 
are divergent in e = {4 — D)/2 at intermediate steps. Obviously, these divergences cancel 
when the intermediate expressions are combined to give flnal amplitudes. An example of 
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this behavior is given by Eqs. ( |C9| ) and ( |(J.l(j| ), where the divergent part of the Bo functions 
cancels among the different contributions to the triangle graphs. 

Given the large number of contributing Feynman graphs, it is most convenient to give 
analytic results for the scattering amplitudes for fixed polarizations of the external quarks and 
gluons. These amplitudes are then evaluated numerically, instead of using trace techniques 
to express polarization averaged squares of amplitudes in terms of relativistic invariants. We 
proceed to derive explicit expressions for these amplitudes. 

3 Notation and matrix elements 

Within the SM, the effective interaction of the Higgs boson with gluons is dominated by 
top-quark loops because the top Yukawa coupling, ht = rrit/v with v = 246.22 GeV, is much 
larger than the Hbb coupling. In the following we only consider top-loop contributions. 
All the H + 4 parton amplitudes, at lowest order, are then proportional to hf g'^, where 
Qs = \/4:7ias is the SU(3) coupling strength. It is convenient to absorb these coupling 
constants into an overall factor 

F = h,-^^4m, = 4^al, (3.1) 

where we have anticipated the loop factor l/lGvr^ and the emergence of an explicit factor, 
Arrit, from all top-quark loops, which results from the compensation of the chirality flip, 
induced by the insertion of a single scalar Hti vertex. 

3.1 qQ — > qQH and qq — > qqH 

The simplest contribution to if -|- 2 jet production is provided by the qQ qQH process 
depicted in Fig. |l|(a). Other four-quark amplitudes are obtained by crossing. We neglect 
external fermion masses and use the formalism and the notation of Ref. [^] for the spinor 
algebra. For a subprocess like 

^i) + Q{Pai h) q{P2i ^2) + <5(p3' ^3) + H{P) , (3.2) 

with 

P?=P2=P3=P4 = 0, P' = m^, (3.3) 
each external ( ant i-) fermion is described by a two-component Weyl-spinor of chirality r = 

Here pj, and ii denote the physical momentum, the helicity and the color index of the quark 
or anti-quark, and the sign factor Si allows for an easy switch between fermions [Si = +1) 
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and anti-fermions {Si = —1). The quark-gluon vertices of Fig. 0(a), including the attached 
gluon propagators, are captured via the effective quark currents 

J^l = xlM)i^')rXaM)j^^^^ = 5..<r, (2| {a^^)r |1) ^^^^^ ^^.S) 

and 

J'a, = 5.4.3 (4| (a^). |3) ■ (3-6) 

Here we have used hehcity conservation via the assignments t = cfi = (J2 and r = cts = cr4, 
respectively, and the sign factors for the fermions provide an easy connection between the 
gluon momenta qi = P2 — Pi and q2 = Pi — Ps going out of the top-quark triangle and the 
physical quark momenta = SiPi. Finally we have used the shorthand notation 



ll) 


= Xai(Pi) , 




(2| 


= XU{P2) , 




|3) 


= XaAPs) , 




(4| 


= xi^pd , 


(3.7) 



and {cr^)± = (l,±cr) is the reduction of Dirac matrixes 7^^ into the two-component Weyl 
basis. Since the quark currents J12 and J43 are conserved, we have 

Ji'2 {Pi - P2), = 0, J^, (P3 - Pa)^ = . (3.8) 

The Weyl spinors and the currents J21 and J43 are easily evaluated numerically 0. The 
scattering amplitude for different flavors on the two quark lines is then given by 

where the overall factor 



F«« = s, S2 Ss 4^p^,pIpIpI F (3.10) 

includes the normalization factors of external quark spinors. The = A° /2 are the color 
generators in the fundamental representation of SU(A^), = 3. As detailed in Appendix |C|, 
the tensor T^y{p, q) can be written as 

r^'^(p,g) = FT{p\q\{p + qf) {p-qg'-'' -p'qn 

-t- g^ (p + g)') [qyg^"" -p^qi^q" -q^p^p" +p-qp^'q'') . (3.11) 

Analytic expressions for the scalar form factors Ft and Fl are given in Eqs. (|C.9|) and ( |C.1CI| ). 

The scattering amplitude for two identical quarks is obtained from the result above by 
including Pauli interference, which results from interchanging quarks 2 and 4, 

AQQ _ AQQ fa j.a _ Aqq j.a /o Tn\ 

— -^2143 '^iaii ''Mis -^4123 '^Mii '-1213 ■ [O.IZ) 



The squared amplitude, summed over initial- and final-particle color, becomes 

E = (I Ai43r + I Amr) + 2 Re (Ams^Iiss) (3-13) 

col 

The squared amplitude for the qQ qQH process can be read from Eq. (|3.13|) by putting 
Al23 = 0. 



3.2 qg qgH 



P2 < Pj 




'I 

i 



P2 e?-) _ 




(gl ^ g2) 



(1) 



(2) 



(3,4) 




a'P2-pi 

P2 <^ Pl 




grP2-Pl 
P2 ^ Pl 



'i' ^ qi+q2 i2 




P2 |P^-P' Pl 



(5) 



(6) 



(7) 



qi 



-H Oliipjlil. 

qi 



q2 q2 




M,qjL(b 
qi 



3' q3=P2-Pl 
P2 ^ Pl 



q2 



q3=P2-Pi 

P2 ^ Pl 



3' q3=P2-pi 
P2 °i Pl 



(8) 



(9) 



(10) 



Figure 2: Feynman graphs contributing to the process qg^qgH. Graphs (3) and 
(4) are the same as (1) and (2), hut with gluon labels interchanged. No distinction 
is made between the two orientations of the fermion arrow on the top- quark loop, 
because they are related by Furry 's theorem. 

Twenty distinct Feynman graphs contribute to the process 

g{qi,ai) + qiPi,ii)^giq2,a2) + g(p2>^2) + H{P) , (3.14) 
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and crossing related processes (the color index of the external gluons is indicated with Oj). 
However, pairs with opposite directions of the top-quark fermion arrow in the loop are related 
by charge conjugation and will not be counted separately in the following (see Appendixes P 
and0). 

The resulting ten Feynman graphs are depicted in Fig. |^. Following Ref. all gluon 
momenta are treated as outgoing. For the specific process of Eq. ( |3.14|) we then set qi = —Tji 
and q2 = +q2, with 

(3.15) 



_2 —2 —2 —2 n 

Pl= P2 = 11 = (l2 = 0, 



Feynman graphs with a triangle insertion in an external gluon line, i.e. with one light-like 
gluon attached to the top-quark triangle, receive contributions from a single form factor. Fx, 
only (see Eq. (|C.12|) ). These simplifications are best captured by replacing the polarization 
with the effective polarization vectors 



vectors e^^ = ef (o. 



e>;H = FT{0,{P + q.)\P' 

^ ^ {P + qi 

for gluons i = 1,2, with Ft given in Eq. ( |C.10|) . 



(gfe,-P-erg,-P) 



(3.16) 



External quark lines are handled as for the qQ^qQH amplitudes. Feynman graphs 5 
through 10 are proportional to the quark current J21 defined in Eq. ( |3.5|) . Spinor normal- 
ization factors are absorbed into the overall factor 



F^^ = -S,S22^TfiTP2F5^,,,. 
Using the shorthand notation 



1 



ij)2 + qiY ' 
1 



(3.17) 

(3.18) 
(3.19) 



for the 2-component Weyl spinors describing emission of a gluon next to an external quark, 
we arrive at a compact notation for the contributions to the qg — qgH scattering amplitude 



a2^ 



(2|(^,^).Jg2,l) + (2,gi|(^2/f)<xJl) 



+ 



(t-t-),^,^ {2\{^,^U\q„l) + {2,q2\ (^1/^)0-1 |1) + 



+ 2(^21 ■ q2 eiH ■ £2 - J21 ■ eiH {P2 - Pi) ■ £2 - J21 ■ £2 92 ■ ein) 
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2(^21 • qi e2H ■ ei - J21 ■ e2H {P2 - Pi) ■ ei - J21 ■ eiQi- e2H , 



X 



ei ■ ^2 (^2 - qi)^ - 2 g2 ■ ei £2/. + 2 gi ■ £2 



R ^/ii M J^J.3 



(3.20) 



where = Ft ((gi + 52)^ , (pi — ^2)^ 7 -P^) , and analogously for Fl. 



The contributions of the box diagrams enter in the last line of Eq. ( p. 20 ) via the tensor 



B 



B 111 ^i2f^3 {91^^2,93} (with ^3 = P2 — Pi)- Gauge invariance and Bose symmetry of 
the gluons limit the relevant structure of these box contributions to just two independent 
scalar functions, as we will now show. 



It is easy to see that the qg qgH amplitude A'^^ is invariant under the replacements 
-^€1 + KiQi and €2 — > €2 + H2q2, for arbitrary constants Kj.Q By proper choice of these 
constants, the polarization vectors of the two on-shell gluons can be made orthogonal to 
both qi and ^2- This can be seen by introducing a convenient basis of Minkowski space, 
composed of the vectors , ^2 and 



x'' = q2- gs qi+qi- gs q2 -qi- q2 gs , 
y" = e^"^^gi,g2/3g3p. 



(3.21) 
(3.22) 



The vector y'^ is orthogonal to all momenta occurring in the boxes (gi,g2 and q^) while 
is orthogonal to gf, g2 and y'^. More precisely 



X ■ X 

y ■ y 

X ■ gi 

y ■ qi 



detQs = gi ■ g2 gl gi ■ g2 - 2 gi ■ gs g2 ■ gs 
detQs , 

detQs 



X ■ q2 

y ■ q2 



, X • gs = 

gi • g2 

0, i/-gs = 0, y-x = 



(3.23) 
(3.24) 

(3.25) 

(3.26) 



where detQs denotes the Gram determinant of the box, i.e. the determinant of the 3x3 
matrix with elements (Qs)jj = —qi ■ qj- Obviously detQs is symmetric under interchange of 
the three gluon labels. The non-symmetric form given in Eq. (|3.23|) pertains to our particular 
situation where only gs may be off-shell {qf = g| = 0). 

^These gauge invariance conditions provided a stringent test of our numerical programs. 
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In the following we take 



^2 



V-detQa 
V-detQs 



(3.27) 
(3.28) 



as the two independent polarization vectors of each of the on-shell gluons. With this choice, 
we eliminate the two Be contributions in the box tensor -B^i/x2M3 (^^e Eq. (p.6| )) since they 
contain a factor gf^ or qi^^ , which vanishes upon contraction with the polarization vectors. 

We can then write the squared element for qg qgH, summed over the polarization 
vectors of the external gluons in the form 



pol 



VdetQg 



l/l^ffpi \ Aqg\2 , \ Aqg\2 , \ Aqg \2 
I •^xx I ~ I "^xy I ~ I "^vx I ~ I "^uy I 



yx\ 



(3.29) 



where the shorthand A^^y = Alf^^^ x^^ y^^, etc. has been used. Expressions for the contracted 
tensor integrals for the boxes that appear in Eq. ( p.20| ) are given in Eqs. (p.l5|) - (p.2CI|) . 

In addition, the color structure of the qg^qgH amplitude is given by (see Eq. (|3.2CI| )) 
SO that the resulting color-summed squared amplitude takes the form 



col 



qg\2 

21 1 



AN 



2 Re 



Al^ {A?' 



21) 



AN 



(3.31) 



3.3 gg^ggH 

For the process 

9{Qi, ai) + 9iQ2^ 0-2) giQs^ ^3) + giq^, 0,4) + H{P) , (3.32) 

we introduce the outgoing momenta qi, so that qi = —qi, q2 = —q2, = +Q3 and ^4 = +^4, 
where 

t=t = t = t = 0, = ml, (3.33) 

and are the color indices in the adjoint representation carried by the gluons. 

Due to the large number of diagrams and the length of the results, we are not going to 
write explicitly the expressions for the amplitude, but we describe in detail the procedure 
we follow. 

We used QGRAF [|I^ to generate the 49 diagrams for this process. As detailed in Sec. 0, 
these diagrams are obtained by the insertion of a triangle, a box or a pentagon loop into the 
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tree-level diagrams for gg gg scattering, so that we can write the un-contracted amplitude 
in the "formal" way 

19 18 12 

{•A^^)ap^y6 = '^ajJ'yS + X! '^i ^ajJ-yS + X! ^alS-yS- (3.34) 

1=1 1=1 i=l 

where the tensor functions T*, and are Feynman diagrams that contain a triangle, a 
box and a pentagon fermionic loop, while cj , cf and cf are the respective color factors. This 
amplitude will be contracted with the external polarization vectors of the gluons, e^, to give 

A^' = 6? 6^ el ei iA^%^^, ^ (A^X.,^^., ■ (3-35) 

We used Maple to trace over the Dirac 7 matrixes and to manipulate the expressions. Since 
the Higgs couples as a scalar to the massive top quark in the loop, the resulting trace in the 
numerator of the generic n-point function has at most n — 1 loop-momentum factors. Using 



the tensor-reduction procedure described by Passarino and Veltman [|1T| , we can express the 
triangle and box one-loop tensor integrals in terms of the external momenta gf and in terms 
of the scalar functions Cij (z = 1 . . . 2, j = 1 . . . 4) and Dij [i = 1 . . .3, j = 1 . . . 13). For 
speed reasons, we preferred not to express directly the final amplitude in terms of scalar 
triangles and boxes (the Co and Dq functions in Passarino- Veltman notation), in our Monte 
Carlo program. Expressions for the amplitude written in terms of the Cq and Dq integrals 
are considerably larger than the result where we keep the Cij and Dij functions. 

In dealing with diagrams with a pentagon loop, we worked directly with the dot products 
in the numerator. In fact, the generic tensor pentagon appearing in gg—^ggH scattering 
has the form 

d^k {fc", k"ki^, k^k^k\ k'^k'^k^k^} 



_ ^2j p _^ _ ^2j p _^ p^^p _ ^2j p _^ p^^^Y _ ^2j p _^ ^1234)^ - ml 

(3.36) 

where pij = Pi+Pj, and similar ones for piji and pijin, and where the set ofpi, {pi,P2,P3,P4}, is 
one of the 24 permutations of the external gluon momenta {qi, q2, gs, q^}. These permutations 
are reduced to 12, once Furry's theorem is taken into account. The generic scalar five-point 
function, that is Eq. (p.36|) with a 1 in the numerator, will be indicated with i?o(pi,P2,P3,P4)- 

The tensor indices appearing in the numerator are always contracted with one of the 
following: 

1. the metric tensor g'^'^ 

2. one of the external momenta gf 
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3. one of the external polarization vector ef . 



The procedure we have used in these three cases is the following. 



1. Every time there is a A;^ = k°'k^gaf3 product in the numerator, we write it as 



1 2 2 

k — 



(3.37) 



and the first term is going to cancel the first propagator, giving rise to a four-point func- 
tion, while the last one will multiply the rest of the tensor structure in the numerator, 
that now has been reduced by two powers of the loop momentum k. 

2. We rewrite every scalar product of the type {k ■ g,) in the numerator as a difference of 
two propagators, using the identity 



k- qi = ]:[[{k+p + Qi) 



m. 



{k + pf 



m. 



2 ■ p} , 



(3.38) 



where p is an arbitrary momentum. The first two terms in the sum are going to cancel 
the two propagators adjacent to the external gluon leg with momentum g^, while the 
last one will contribute to the rest of the tensor structure in the numerator, now reduced 
by one power of k. 

3. When the dot product {k ■ €{) appears in the numerator, we choose the four external 
gluon momenta as our basis of Minkowski space. Hence we can expand the external 
polarization vector Cj as 



4 

E 



(3.39) 



where the coefficients ej,- are computed by inverting the system of equations 



ei-qk = (24),fc, hk 



where the elements of the matrix Q4 are given by 



(3.40) 



(3.41) 



Using Eq. ( p.39| ), we can rewrite every scalar product of the form (k-ei) in the numerator 
of the Feynman diagrams as a sum over {k ■ qj), that we handle in the same way as 
shown in Eq. ( p. 381) . 

The matrix Q4 is singular when the four gluons become planar, i.e. when the four gluon 
momenta cease to be linearly independent. This can easily be seen in the center-of-mass 
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frame of partons 1 and 2, where we can write 



q^ = E(1,0,0,1), 

q2 = E (1,0, 0,-1), 

g3 = (1, sin 6*, 0, cos 6*) , 

g4 = £'4 (1, sin cos 0, sin6'' sin 0, cos6'') , 

and the determinant of the matrix Q4 becomes 

det Q4 = -4.E^ElEl (1 - cos^ 6') (l - cos^ ^) (l - cos^ 0) . (3.42) 

The cuts imposed on the final-state partons (see Eq. ( ^.l]) ) avoid the singular region 
when one of the final-state partons is collinear with the initial-beam direction (singu- 
larities in 9 and 6'). 

The singularity in cf) is un-physical, and the final amplitude should be finite near 
the singular points (ps = 0, vr. In our FORTRAN program, when the Monte Carlo 
integration approaches the singular points in 0, we interpolate the value we need from 
the values of the amplitude in the points = 0s =t 0.01 vr. We have checked that the 
interpolated amplitude differs from the exact value of the amplitude by less than 1%, 
in the non-singular region. 

By iterating this reduction procedure, we can write the contracted amplitude for gg ggH 
scattering in terms of 

- twelve EQ{pa,Pb,Pc,Pd) functions, i.e. the scalar five-point functions computed in dif- 
ferent kinematics; 

- twelve permutations of Dij functions with argument Dij{pa,Pb,Pc), with a < c, six 
Dij{pa,Pb +Pc,Pd) with a < d and twelve Dij{pa,Pb,Pc +Pd), together with the corre- 
sponding Do functions; 

- three Cij{pa+Pb,Pc+Pd), and four Cij{pa,Pb+Pc+Pd), together with the corresponding 
Co functions. 

As usual, the set {Pa,Pb,Pc,Pd} is chosen in the group of permutations of the gluon momenta 
{gi, ^2, ^4}- To reduce the number of the Dij and Cij functions to this independent set, 
we made use of some identities between these functions, that we collect in Appendix |B[ 

As stated in Eq. ( p.34| ), we can study the color factors of the gg^ggH process, by 
dividing the full amplitude into three different classes, according to the number of fermionic 
propagators in the loop. We first discuss the color factors of the diagrams containing a 
pentagon loop. 
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Diagrams with a pentagon loop 

The contribution from the sum of charge-conjugated pentagon diagrams is proportional 
to the sum of two color traces with four t matrixes (see Eq. ([E.3|) ). From the invariance 
property of the trace under cyclic permutations, we have only (4 — 1)! = 6 independent 
traces, from the permutation of the four gluon indices. These six color traces combine 
together as in Eq. ( Pi).3|) to give rise to three independent color structures 

C2 = Ti {en^'H^H''^) + Tt {t^H''H''H''^) (3.43) 
C3 = Tr {t''H''H''H'''') + Tr {t^H^H^H"^) . 

The Cj color coefficients are real. In fact, using the identity (the sum over the repeated 
index is understood) 

1 1 
AN 8 ^ ^ ^ ^ 

where / is the (totally antisymmetric) SU(A^) structure constant and d is the totally 
symmetric symbol, we can write, for example. 



1 



2 

N 



(3.45) 



and similar ones for C2 and C3. Since / and d are real constants, the c, are real too. 
A few useful identities can be derived if we take the differences of the q 

Cl - C2 = -^/"i^^'r'""' =^ = 2(C2 -Cl), 

_ = _\pia,l p2a;l ^ fa,a,lj^a2a,l ^2{ci-C3), (3.46) 

C2 - C3 = -^/''^''^V*'""' =^ Z''^"^'/'''"'"' = 2 (C3 - C2). 



Note that the differences of the Cj in the system (|3.46| ) automatically embodies the 



Jacobi identity: by summing the three expressions, we have 

j}aia2l j^a^ail _|_ j:aia4l j!a2a3l _|_ jaia-^l ^0402? g 4'^^ 

Diagrams with a box loop 

These diagrams all contain a three-gluon vertex together with the quark loop. Since 
the sum of the charge-conjugated boxes is proportional to the structure constant / 
(see Eq. ( |D.3|) ), the final color factors accompanying these diagrams are a product of 

two /'S, such as f^i^^l faaa^l _ 



With the help of Eq. ( p.46| ), we can express these products in terms of differences of 
the c,- color factors. 
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Diagrams with a triangle loop 

The same argument can be used to show that the color structure of all the diagrams 
with a three-point function insertion are proportional to the product of two structure 
constants, that are then converted to differences of Cj color factors, using the identities 
of Eq. ( CT) . 



Since all the color structures of the diagrams contributing to the gg—^ggH process can 
be written in terms of the q color structures of Eq. ( p.43| ), we can then decompose the full 
amplitudes in the following way 

3 

j,99 = Y,CiAf. (3.48) 

1=1 

The sum over the external colored gluons of the squared amplitude becomes 

col *ii=l col 



where we have taken into account the fact that the q are real (see Eq. ( p.45| )). Using 
Eq. (|3^) , one finds 

V- (A^2 _ _ 2A^2 ^ 6) 
Ci = 2^CiCi = , no summation over i (3.50) 

col 

r - \- (iV^ - 1) [N^ - 3) . 

C2 = }_^CiCj = , (3.51) 

col 

and we finally get 

col i=l 



4 Checks 

We were able to perform two different kinds of checks on the analytic amplitudes we 
computed: a gauge-invariance and a large-m^ limit check. 



4.1 Gauge invariance 

Gauge invariance demands that the amplitudes should be invariant under the replacement 
— > Cj + Kiq-i, for arbitrary values of Kj. This implies that for the qg — >■ qgH process we must 
have (see Eq. (Q)) 



when qi-ei = 0, i = 1...2, (4.1) 
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and for gg ^ ggH (see Eq. ( |3.35D ) 










when qi ■ ei = 0, 



1...4 



(4.2) 



We checked the gauge invariance in Eqs. (4J.) and (4^) both analytically and numerically 
(in the final Fortran program). 



Using Eq. (|3.2CI|) , it is straightforward to check that the system ( [4.1|) is satisfied. To 
check gauge invariance for the gg ggH process, we wrote the contracted amplitude of 
Eq. ( |3.35|) in terms of scalar pentagon (£^0)5 box (-Do), triangle (Co) and bubble (i?o) integrals, 
keeping the space-time dimension D arbitrary. This means that we expressed the Cij and 
Dij functions in terms of Bq, Cq and Dq "master" integrals. The coefficients of these scalar 



integrals are then functions of scalar products (qi-qj), {qi-e 
introduced in Eq. (|3.39| ). 



and of the e,i coefficients 



Since the two-point functions Bq are divergent in e = (-D — 4)/2, and since the total 
amplitude must be finite, these poles must cancel. In fact, the factors multiplying the Bq 
contributions are proportional to (1^ — 4), so that only the pole coefficient of Bq contributes 
to the finite amplitude (see comment after Eq. ( |C10D ). 

To implement the gauge invariance check with respect to the polarization vector e^, as 
described in the system (^.2|), we make the replacement 



0, jV^, 



(4.3) 



in Eq. ( |3.35 ), and we impose the orthogonality condition ek ■ qu = 0, that constrains the euj 
coefficients to satisfy the identity (see Eq. (|3.39|) ) 



4 



■ W 



0, 



(4.4) 



We have checked gauge invariance in two different ways. 



1. Suppose that instead of considering the QCD process gg-^ ggH, we consider the QED 
analogue, 77 '-yjH. In this scenario, all the diagrams with a three- or a four-gluon 
vertex are no longer present: the only surviving diagrams are the ones containing a 
pentagon loop, with no color structure associated. The amplitude, not contracted with 
any external photon polarization vectors, is (see Eq. ( |3.34|) for comparison) 

2 = 1 



16 



The gauge invariance of this expression allows us to check the correctness of the tensor 
reduction of the pentagon diagrams only. We have contracted Eq. ([4.5|) with the 
polarization vectors of the photons and we have applied the tensor reduction procedure 
previously described. Instead of expressing the results in terms of Cij and Dij functions, 
we have expressed these coefficient functions in terms of Sq, Cq, Dq and Eq scalar 
integrals, keeping the space-time dimension D arbitrary. Since these scalar integrals 
form a set of independent functions, we expect the coefficients of these integrals to be 
zero, in order to fulfill the gauge- invariance test. Note that we have considered the 
twelve scalar five-point functions Eq as independent from the four-point functions Dq, 
that is we have not used Eq. ( |A.1| ). This is indeed the case, since, in arbitrary D 



dimensions, the scalar pentagon cannot be expressed as a combination of scalar boxes 
only, so that it is really an independent integral. 

Finally, we have checked that our full QCD amplitude satisfies the four identities in 
the system ( [4.2|) . Since the amplitude can be split into three different contributions 
according to the three independent color factors Cj (see Eq. ( p.4S )), this means that 



not only the full amplitude is gauge invariant, but that the three sub- amplitudes Af^ 
are separately gauge invariant, and satisfy a system of equations similar to 



4.2 Large-mi limit 

The amplitudes for Higgs plus two partons agree in the large-m^ limit with the corresponding 



amplitudes obtained from the heavy-top effective Lagrangian |]12|. This check was done 
numerically by setting = 3 TeV. We found good agreement with the — ^ oo results, 
within the statistical errors of the Monte Carlo program, for Higgs boson masses in the 
range 100 GeV < mn < 700 GeV. 



5 Applications to LHC physics 

The gluon-fusion processes at O^a^), together with weak-boson fusion {qq^qqH pro- 
duction via t-channel exchange of a W or Z), are expected to be the dominant sources of 
+ 2 jet events at the LHC. The impact of the former on LHC Higgs phenomenology is 
determined by the relative size of these two contributions. However, the gluon-fusion cross 
sections ioi H + 2 jet events diverges as the final-state partons become coUinear with one 
another or with the incident beam directions, or as final-state gluons become soft. A minimal 
set of cuts on the final-state partons, which anticipates LHC detector capabilities and jet 
finding algorithms, is required to define an + 2 jet cross section. Our minimal set of cuts 
is 

PTj > 20 GeV, \rij\ < 5, % > 0.6, (5.1) 
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where pTj is the transverse momentum of a final state parton and Rjj describes the separation 
of the two partons in the pseudo-rapidity rj versus azimuthal angle plane 



(5.2) 
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Figure 3: H + 2 jet cross sections in pp collisions at ^/s = 14 TeV as a function of 
the Higgs boson mass. Results are shown for gluon-fusion processes induced by a top- 
quark loop with mt = 175 GeV and in the — > oo limit, computed using the heavy- 
top effective Lagrangian, and for weak-boson fusion. The two panels correspond to 
two sets of jet cuts: (a) inclusive selection (see Eq. ( [5. l\ ) ) and (b) WBF selection 



(Eqs. and (j^)). 



Expected + 2 jet cross sections at the LHC are shown in Fig. ^ as a function of the 
Higgs boson mass, irtH- The three curves compare results for the expected SM gluon-fusion 
cross section at mt = 175 GeV (solid line) with the large-mt limit (dotted line), and with 
the WBF cross section (dashed line). Error bars indicate the statistical errors of the Monte 
Carlo integration. Cross sections correspond to the sum over all Higgs decay modes: finite 
Higgs width effects are included. 



In all our simulations, we used the CTEQ4L set for parton-distribution functions ||13 . 
Unless specified otherwise, the factorization scale was set to /i/ = y/pTiPT2- Since this 
calculation is a LO one, we employ one-loop running of the strong coupling constant. In 
Fig. I we fix as = as{Mz) = 0.12. 



The left panel in Fig. ^ shows cross sections within the minimal cuts of Eq. ( p.l|) . The 
gluon-fusion contribution dominates because the cuts retain events with jets in the central 
region, with relatively small dijet invariant mass. In order to assess background levels for 
WBF events, it is more appropriate to consider typical tagging jet selections employed for 
WBF studies [^. This is done in Fig. ^ (b) where, in addition to the cuts of Eq. (pTT|), we 
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Figure 4: H + 2 jet contributions to the cross section in pp collisions at ^/s = 
14 TeV as a function of the Higgs boson mass. Results are shown for the different 
contributions to the gluon-fusion process (gg,qg and qq amplitudes) using (a) the 
inclusive cuts of Eq. l \5.J\ ) and (b) the WBF cuts of Eqs. ( |5.i| j and ( [5. 4 )- 



require 

kji - rij2\ > 4.2, rjji ■ r]j2 < 0, rnjj > 600 GeV, (5.3) 

i.e. the two tagging jets must be well separated, they must reside in opposite detector hemi- 
spheres and they must possess a large dijet invariant mass. With these selection cuts the 
weak-boson fusion processes dominate over gluon fusion by about 3/1 for Higgs boson masses 
in the 100 to 200 GeV range. This means that a relatively clean separation of weak-boson 
fusion and gluon-fusion processes will be possible at the LHC, in particular when extra 
central-jet- veto techniques are employed to further suppress semi-soft gluon radiation in 
QCD backgrounds. We expect that a central-jet veto will further suppress gluon fusion with 
respect to WBF by an additional factor of three M . 



A conspicuous feature of the if +2 jet gluon-fusion cross sections in Fig. | is the threshold 
enhancement at m/f ~ 2m(, an effect which is familiar from the inclusive gluon-fusion cross 
section. Near this "threshold peak" the gluon-fusion cross section rises to equal the WBF 
cross section, even with the selection cuts of Eq. ([5.31). Well below this region, the large 
rrit limit provides an excellent approximation to the total if + 2 jet rate from gluon fusion, 
at least when considering the total Higgs production rate only. Near top-pair threshold the 
large rrit limit underestimates the rate by about a factor of 2. 

A somewhat surprising feature of Fig. ^ (b) is the excellent approximation provided by 
the large rrit limit at Higgs boson masses below about 200 GeV. Naively one might expect the 
large dijet invariant mass, rrijj > 600 GeV, and the concomitant large parton center-of-mass 
energy to spoil the rrit ^oo approximation. This is not the case, however. As shown in 
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Figure 5: Renormalization- scale dependence of the total cross section for H plus 



two jet production with the inclusive cuts of Eq. (5.1). The renormalization scale 
fir = (.fJ-o is varied in the range 1/5 < ^ < 5. The five curves correspond, from top to 
bottom, to the following choice of fiQ.- the geometric mean of the transverse momenta 
of the two jets, the Z mass, the invariant mass of the two jets, the geometric mean of 
the two invariant masses of the Higgs and the jets, and the partonic center- of-mass 
energy. 



Ref. [§], the large rrit limit works well in the intermediate Higgs mass range, as long as jet 
transverse momenta stay small: pxj ^ rnt. 

In Fig. ^ we have plotted the individual contributions to the gluon-fusion cross section 
which are coming from the gg ggH, qg — ^ qgH and qq — * qqH sub-processes, including all 
crossed processes for each of the three subgroups. Results are shown after imposing the 
inclusive cuts of Eq. (|5.1|) (left panel) and the WBF cuts of Eqs. (|5.1| ) and (|5.3|) (right 
panel), with rrit = 175 GeV, so that the sum of the three curves in each panel add up to the 
solid line curve in Fig. ^ External gluons dominate in the inclusive-cut case (left panel): 
final-state gluons tend to be soft and initial gluons preferably lead to soft events due to the 
rapid fall-off of the gluon parton distribution function, g{x, fif), with increasing x. When the 
rrijj > 600 GeV constraint of the WBF cuts is imposed, the gluon contribution dies rapidly, 
as is shown in Fig ^ (b). 

The results shown in Fig. ^ raise two questions, which we intend to answer in the fol- 
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Figure 6: Factorization- scale dependence of the total cross section for H plus two 
jet production with the inclusive cuts of Eq. ( ^5. The factorization scale fJ.f = ^/Uq 
is varied in the range 1/5 < ^ < 5. The three curves correspond to the following 
choice of hq: the geometric average between the transverse momenta of the two jets, 
the invariant mass of the two jets and the partonic center-of-mass energy. 

lowing: i) what are the uncertainties in the prediction of the + 2 jet cross section, and 
ii) which distributions are the most effective in distinguishing gluon-fusion and weak-boson 
fusion contributions to H + 2 jet events? 

In order to asses the sensitivity of the gluon-fusion cross section to higher order QCD cor- 
rections, we have plotted in Fig. |^ the total cross section for several choices of the renor- 
malization scale (the factorization scale has been kept at /i/ = ^"PtiVti)- We have fixed 
A|jg = 254 MeV, so that, as{Mz) = 0.12, with nj = 5 active flavors. We have chosen five 
different scales: /xq = ^JPTlPT2■, Mz, rrijj, y/rriHji mHj2 and V^, i.e. the geometric average 
between the transverse momenta of the two jets, the Z mass, the invariant mass of the 
two jets, the geometric mean of the two invariant masses of the Higgs and the jets and the 
partonic center-of-mass energy. For every event generated by our Monte Carlo, we have 
computed the running of the coupling constant agifJ^r) at the values fir = ^fio, where ^ was 
allowed to vary from 1/5 to 5. We can see that the renormalization-scale dependence is very 
strong, mainly due to the fact that this is a leading-order calculation, at order a^. 

What is the "most natural" scale for as is an unresolved issue. The good agreement 
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between the complete result and the nit ^oo one (away from threshold) implies that the 
cross section is dominated by Feynman diagrams with a gluon exchange in the t channel. 
These diagrams contain a triangle loop that couples the t-channel gluon with the Higgs. For 
this reason, it seems reasonable to make the replacement 



otsiPTi) as{PT2) alimu)- 



(5.4) 



With this choice for the strong coupling constant we have a total cross section (within the 
cuts of Eq. (p.lD) of about a = 9.6 pb, which sits in between the two values computed 
using ag{Mz) and af(mjj) (see Fig. |]). Dismissing the extreme choice fXr = (which is ill 
defined at higher orders), and allowing for the conventional factor of 2 variation of ^, Fig. ^ 
suggests an uncertainty of the gluon fusion if + 2 jet cross section of about a factor 2.5 as 
compared to the central value of 9.6 pb obtained with the central choice of Eq. ( ^.4|) . 



Keeping fixed the renormalization scale at Mz, we have collected in Fig. |^ the results for 
the factorization-scale dependence of the total cross section. The factorization scale /i/ = ^/iq 
was allowed to vary in the range described by 1/5 < ^ < 5, where /xq was taken equal to 
the geometric average of the transverse momenta of the two jets, the invariant mass of the 
two jets and the partonic center-of-mass energy. Compared with the variation with respect 
to the renormalization scale, we see that the dependence on the factorization scale is almost 
negligible: in Fig. 13 the if + 2 jet cross section varies between 9.2 pb and 10.2 pb. 



In the following we use the renormalization scale choice of Eq. ( ^.4|) and set /i/ 
y/PTi Pt2- We take rriH = 120 GeV throughout, as a characteristic Higgs boson mass. 
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Figure 7: Dijet invariant-mass distribution of the two final jets for gluon-fusion 
(solid) and WBF (dashes) processes. Left panel (a): inclusive cuts of Eq. ( [5. right 
panel (b) WBF cuts of Eqs. ( [5. J| j and l{5.S^ , where we have suppressed the constraint 



nijj > 600 GeV. 
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Figure 8: Rapidity separation of the two final jets for gluon-fusion (solid) and WBF 
(dashes) processes. Left panel (a): inclusive cuts of Eq. ( [5. right panel (b) WBF 
cuts of Eqs. ( |5../| j and ^5.3( j, where we have suppressed the constraint \r]ji—rjj2\ > 4.2. 



Turning now to the issue of differentiating between gluon fusion and WBF processes, 
the prominent characteristics to be considered here are the jet properties. Figure |^ shows 
the dijet-mass distribution for gluon-fusion and WBF processes, using the inclusive and 
the WBF cuts. In this last case, we have suppressed the constraint nijj > 600 GeV, in 
order to access the region of small dijet invariant mass. In both panels, the high dijet mass 
region [rrijj > 1 TeV) is dominated by WBF. The significantly softer dijet-mass spectrum 
of the gluon-fusion processes is characteristic of QCD processes, which are dominated by 
external gluons, as compared to quarks in WBF processes (see Fig. ^ and comments about 
it). Figure |^ also shows the rrit ^oo dijet-mass distributions (dotted curves), that are 
almost indistinguishable from the rrit = 175 GeV result: large dijet invariant masses do not 
invalidate the ^ oo limit as long as the Higgs boson mass and the jet transverse momenta 
are small enough, less than mt in practice. 

A characteristic of WBF events is the large rapidity separation of the two tagging jets, 
a feature which is not shared by if + 2 jet events arising from gluon fusion. The rapidity 
separation of the jets is shown in Fig. ^, for both gluon-fusion (solid) and WBF (dashes) 
processes. The two panels correspond to the inclusive cuts of Eq. (|5.1|) and to the stricter 
WBF cuts of Eq. (|5.3|) , where we have suppressed the constraint \riji — 77^2! > 4.2, in order to 
have access to the entire Arjjj range. The jet separation cut, \riji — r^^^l > 4.2, is one of the 
most effective means of enhancing WBF processes with respect to gluon fusion. The small 
dip in the gluon-fusion distribution at small Arjjj is a consequence of the cut Rjj > 0.6. 

A second jet-angular correlation, which allows to distinguish gluon fusion from weak- 
boson fusion, is the azimuthal angle between the two jets, (pjj. The distributions for gluon- 
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Figure 9: Azimuthal-angle distribution between the two final jets, with the WBF 
cuts of Eqs. ( [5. 1\) and (\5.3J. Results are shown for gluon-fusion processes induced 
by a top- quark loop with rrit = 175 GeV and in the mt^co limit, computed using 
the heavy-top effective Lagrangian, and for weak-boson fusion. 



fusion and WBF processes are shown in Fig. In the WBF process qQ — > qQH, the matrix 
element squared is proportional to 

I^WBFp oc — ^ — ^sm^,-, (5.5) 

and is dominated by the contribution in the forward region, where the dot products in 
the denominator are small. Since the dependence of m|j on (pjj is mild, we have the flat 
behavior depicted in Fig. |^. The azimuthal-angle distribution of the gluon-fusion process 
is instead characteristic of the CP-even operator HG^i^G^'^, where G^j,^ is the gluon field 
strength tensor [|I5|. This effective coupling can be taken as a good approximation for the 
ggH coupling in the high-m^ limit. Note that the large-mj limit (dotted line) is almost 
indistinguishable from the = 175 GeV result (solid line). 

Finally, in Fig. we show the transverse-momentum distribution of the Higgs boson in 
gluon-fusion (solid lines) and in WBF (dashes lines) processes, with the inclusive selection 
of Eq. i ^A\} . Within these cuts, both differential cross sections peak around a value of 
pj^ff ^ 50 GeV. Note, however, that, while the peak position of the WBF distribution is 
largely tied to the mass of the exchanged intermediate weak bosons, the peaking of the gluon 
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fusion processes occurs just above 40 GeV, which is a direct consequence of the pxj > 20 GeV 
cut of Eq. (13). 




6 Conclusions 



In the previous sections, we have provided the results of the calculation of H + 2 jet 

cross section, including the full top-mass dependence. For the quark-quark and quark-gluon 
scattering amplitudes we have found very compact analytic expressions. The gg^ggH 
amplitudes, which include pentagon loops, are more complex and available analytically, as 
Maple output, and numerically, in the form of a FORTRAN program. 

Numerical investigations of the resulting cross sections at the LHC provide many inter- 
esting insights. With minimal jet-selection cuts (parton separation of Rjj > 0.6 and jet 
transverse momenta in excess of 20 GeV) the gluon-fusion induced cross section is sizable, 
of order 10 pb for uih = 120 GeV, which corresponds to about 30% of the inclusive Higgs 
production rate. Since our calculation gives the if + 2 jet rate at leading order, it exhibits 
the large renormalization-scale dependence to be expected of an order process. 

As expected, the large-m^ limit provides an excellent approximation to the full de- 
pendence when the Higgs mass is small compared to the top-pair threshold. The Isnge-rrit 
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limit is found to break down for niH > rrit and when jet transverse momenta become large 
{pTj ~ TTT't)- However, large dijet invariant masses do not invalidate the rrit -^oo limit, as 
long as the Higgs boson mass and the jet transverse momenta are small enough, less than 
the top-quark mass in practice. This observation opens the possibility of NLO corrections 
io H + 2 jet production from gluon fusion. Performing the calculation in the large-m^ limit 
would correspond to a 1-loop calculation of a 2 — 3 process. Such a calculation might be 
desirable to reduce systematics errors in the extraction of HZZ and HWW couplings from 
the competing weak-boson fusion processes at the LHC. 

Consideration of the gluon fusion H + 2 ]et rate as a background to WBF studies con- 
stitutes an important application of our calculation. While the overall if + 2 jet rate is 
dominated by gluon fusion at the LHC, kinematic properties of the two processes are suffi- 
ciently different to allow an efficient separation. Gluon-fusion events tend to be soft, with a 
relatively small separation of the two jets. In contrast, the two tagging jets of weak-boson 
fusion events have very large dijet invariant mass, and are far separated in rapidity. Using 
rapidity and invariant-mass cuts, the gluon-fusion cross section can be suppressed well be- 
low the WBF rate. In addition, the azimuthal angle between the two jets shows a dip at 90 
degrees which is characteristic for loop-induced Higgs couplings to gauge bosons . Based 



on our calculation we conclude that a relatively clean separation of weak-boson fusion and 
gluon fusion Higgs plus two-jet events will be possible at the LHC. 
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A Scalar integrals: Cq, Dq and Eq functions 

All the scalar integrals needed for the calculation are finite in D = A dimensions, due 
to the presence of the top-quark mass. No further regulator is required. Scalar triangles 
(Co) and boxes (-Do) have been known for a long time in the literature and efficient 
computational procedures are available [O . Following the procedure outlined in Refs. ITB 



we can express all scalar pentagons as linear combinations of scalar boxes 

5 
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where 



^0 = Do{p2,P3,P4) 

Dl = Do{p4,P3,Pl+P2) 

Dl = Do{p,,p2+P3,Pi) (A.2) 

Dq = Do{pi,P2,P3+Pa) 

Dl = Do{pi,P2,P3) , 



and the matrix T = C ^, with 



Ci, = {r,-r,f -2m\, (A.3) 

and 

ri = 
^2 = Pi 

rs = Pi2 =Pi+P2 (A.4) 

^4 = Pl23 =Pl+P2+P3 

r5 = P1234 =Pl+P2+P3+P4 ■ 

B Relations among Cij and Dij functions 

In this section, we collect a few identities between the Cij and Dij functions. Their 



definition can be found in Ref. |TT[. Please note that we use a (+,—,—, —) metric tensor, so 



that Passarino-Veltman recurrence relations are the same as ours if we make the substitution 
p ■ — P ■ q and — gp^y in their formulae. 



Starting with a three-point vertex function 

d^k f{k) 



m^/'^ [k? — mf] [{k + pY — mf\ [{k + p + qY — mf\ ' 



fB.ll 



where f{k) is an arbitrary function, that, for our purpose, will take the values f{k) = 
1, k", k'^k^, and imposing the equality between this integral and the same integral where 
the integration variable has been shifted according to fc— > — k — p — q, we have 

d^k f{k) 

[{k + p + qy — mt] 

f( — h — n — 

(B.2) 



ZTT^/^ — 777,2] ^(^J^ _|_ p)2 _ j^2j |^^^ _|_ ^ _|_ _ ^2 

d^k fi~k — p — q) 



It^DI2 [^(^/^ _|_ p _|_ _ ^2j \^{^ _|_ ^2 _ ^2j \^2 _ ^2j 
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If f{k) — 1, this identity gives 

Co(p,?) = Co(g,p), (B.3) 

while if f{k) — A;", it gives 

Cii{p,q) = -Ci2{q,p) -Co{q,p), 

Ci2{p,q) = -Ci,{q,p)-Co{q,p). (B.4) 
In a similar way, if f{k) = k°'k^ we get 

C2i(p,g) = C22{q,p) + 2Ci2{q,p) + Co{q,p), 
C22{p,q) = C2i{q,p) + 2Cii{q,p) + Co{q,p), 

C23(p,g) = C'23(g,p) + C'i2(g,p) + Cii(g,p) + Co(g,p), 

C24(p,g) = C24(g,p). (B.5) 

Starting with a four-point function, we derive, in the same fashion. 



^0 


ip,q, 


I) 


= Do{l,q,p), 


Dn 


(p,q, 


I) 


= -Di3(l,q,p) - Do(l,q,p), 


Dl2 


{p,q, 


I) 


= -Di2{l,q,p) - Do{l,q,p), 




{p,q, 


I) 


= -Dn{l,q,p) - Do{l,q,p), 


D21 


(p,q, 


I) 


= D23{l,q,p) + '2Di3{l,q,p) + Do(l,q,p), 


D22 


{p,q, 





^ D22{l,q,p) + 2Di2{l,q,p) + Do{l,q,p), 


D23 


{p,q, 





= D2i{l,q,p) + 2Dn{l,q,p) + Do{l,q,p), 


D24 


(p,q, 





= D26(l, q,p) + Di3(l, q,p) + Di2(l, q,p) + Do(l, q,p), 


D25 


(p,q, 





= D25{1, q,p) + Di3{l, q,p) + Du{l, q,p) + Do{l, q,p), 


D26 


{p,q, 





= D2i{l,q,p) + Di2{l,q,p) + Dii{l,q,p) + Do{l,q,p), 


D27 


{p,q, 





= D27{l,q,p), 


Dzi 


{p,q, 





= -SDi3{l,q,p) -3D23(l,q,p) - 033(1, q,p) - Do(l,q,p), 


D32 


{p,q, 





= -3Di2{l,q,p) -3D22il,q,p) - 032(1, q,p) - Do(l,q,p), 


D3S 


{p,q, 





= -3Du(l,q,p) -3D2i(l,q,p) - D3i(l,q,p) - Do(l,q,p), 


Du 


(p,q, 





= -2Di3(l,q,p) - Di2(l,q,p) -2D2&(l,q,p) - D3g(l,q,p) 








-D23(l,q,p) - Do(l,q,p), 


D35 


{p, q, I) 


= -2 L>i3(/, - D^^(l, q,p)-2 D2^(l, q,p) - 037(1, q,p) 
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-D23{ 




Do{l,q,p), 










i,q,p) - 


2D^2{l,q,p) 


-2D2S,q,p) 


- D3s(l,q,p) 






i,(i,p) - 


Do{l,q,p), 










i,(i,p) - 


2Di,il,q,p) 


-2D25(l,q,p) 


- D35(l,q,p) 






i,(i,p) - 


Doil,q,p), 






Dssip.qJ) = 


-2D,2{l,q,p) 


- Dnil,q,p) 


-2D2,(l,q,p) 


- D3Q(l,q,p) 




-D22{ 


i,q,p) - 


Doil,q,p), 






Dsgip.qJ) = 


-D,2{ 


i,q,p) - 


2D^,{l,q,p) 


-2D2,(l,q,p) 


- D3i(l,q,p) 




-D2l{ 


i,(i,p) - 


Do{l,q,p), 






D3io{p,qJ) = 




i,(i,p) - 


Di2il,q,p) - 


Dii(l,q,p) - D2&(l,q,p) 




—D310 


{i,q,p) - 


- D25il,q,p) - 


- D2i(l,q,p) - 


Do(l,q,p), 


D3iiip,qJ) = 


-D27i 


i,q,p) - 


0313(1, q,p), 






D3i2ip,qJ) = 


-D27{ 


i,q,p) - 


0312(1, q,p), 






^3i3(p,g,0 = 


-D27{ 


i,(i,p) - 


D3ii(l,q,p)- 







(B.6) 



C Tensor and color structure of triangles 
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Figure 11: Two three-point functions connected by charge conjugation. 



The two generic three-point functions depicted in Fig. |TT] have the following expressions 



Tr\qi,q2) 

Tr\qi,q2) 



Arrit J m"^ 
1 r d^k 



Tr 



-7' 



Ml. 



Arrii J in 



-Tr 



- rrit 
1 



-7' 



M2. 



+ ii-mt 
1 



7 



M2. 



1 



-7 



Ml 



1 



(C.l) 
(C.2) 



TTLt if + i^-nit If + ii + i2 - ^t, 

where qi and q2 are outgoing momenta and where we put an overall factor l/(4m() in front 
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to delete a corresponding term coming from the trace over the top quark. Using the charge- 
conjugation matrix C 



(C.3) 
(C.4) 



we can derive (Furry's theorem) 
In addition, 

9i'^mM2(?i>92) = T^iM2(?i>92) = (C.5) 

expresses the gauge invariance of the triangle graphs. The generic tensor structure satisfying 
Eq. (|C3|) is then 

T'-'"' {qi, q2) = Ft (g?, ql (gi + q^f) T^'''{qi, ^2) + {ql ql {qi + 92)') '^^(g^, q,) 
= FtT^'^^ + FlT^^^\ (C.6) 

where (dropping the dependence on external momenta qi and g2, for ease of notation) 

(C.7) 



rplJ.ltM2 



qi-q2 9^'^'-q'^'q!^' , 
ql ql 9'^'^ - ql q^' q^' - ql q'l' + 9i " ^2 q'l'q^' 



(CI 



and 
FLiqlqlQ'] 



2detQ2 



+ 



3g?g2-g 

detQa 

3 gj gi ■ Q 

detQs 



Bo (gi) - Bo (Q) 



Bo (g2) - Bo (Q) 



2 , ^2 3g2g2g2 



4m,' + + ga' + g 



detQs 



Co(gi,g2) + i? , (C.9) 



FriqlqlQ'' 



2detQ2 



Bo (gi) + Bo (ga) - 2Bo (Q) - 2 gi ■ Co (gi, gs) 



+ (gi - g2 [Bo (gi) - (g2) - gi ■ gs i^L • 



(C.IO) 



Here, the negative of Q = gi + g2 denotes the four- momentum of the Higgs boson, detQ2 = 
gigi ~ (gi ■ g2)^ is the Gram determinant, and the terms proportional to i? = —2 are 
pole residues in D = 4 dimensions, originating from contributions proportional to {D — 
4) C24(gi, g2) in the tensor reduction procedure. Please note that even though the Bo func- 
tions are divergent in e = {A — D)/2, the form factors and Ft are finite. 

If one of the external momenta, for example gi, is light-like (real photon or gluon), with 
polarization vector ei, then 



lfJ-2 ^Ml 



gr (gi-g2g2^^-g2'gf ) if gi^ = o 



(C.ll) 
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and as a consequence of the orthogonality gi ■ ei = we have 



T-i/ii/i2 n 



(C.12) 



i.e. the form factor Fl does not contribute when an on-shell gluon or photon is attached to 
the triangle graph. 

The color structure of the sum of the two diagrams in Fig. |lT] is straightforward. In fact, 
both the diagrams have the same color structure and we can write (see Eq. (|C.4| )) 

Tr {fH"^) Tf^"^ (gi, ga) + Tr (t^^t-^i) (g^, q^) = 5"'"^ T^'^"^ (gi, q^) . (C.13) 



D Tensor and color structure of boxes 
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Figure 12: Two four-point functions connected by charge conjugation. 



The two generic four-point functions connected by charge conjugation, depicted in Fig. |12 
have the following expressions 



Bi (gi,g2,g3) 



1 r d^k 



Tr 



Ml. 



1 



-r 



Arrit J in^ \}^ — mt Iji + ^^ — mt 

1 „- 1 



X 



+ 



'12 



Anif J in 



Trl 



-7 



M3. 



7' 

1 



+ 



-7 



M2 



+ 



-7 



Ml. 



'23 



"^t ^ + ^123 ~ "^t 



(D.l) 



where gi, g2 and gs are the outgoing momenta, qij = qi + qj and qijk = qi + qj + qk- The overall 
factor 1 / (4mf ) cancels a corresponding term coming from the trace over the top quark. From 
charge conjugation we have 



Bi (gi,g2,g3) = (gi,g2,g3) = 5 (gi,g2,g3)- 



(D.2) 
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The color structure of the sum of the two diagrams depicted in Fig. ^ is 

Tr {eH^H^^)B^,''"'\qr, q^, gs) + Tr {rH^H^^)B^,"'"'\q,, q^, gs) 

= [Tr(t'^it'^^r3) -Tr(r3t'^^t"i)]B^^'^^^^(gi,g2,g3) = ^^"^"^^'^''''^'(gi, ^2, ga), 

(D.3) 

where we have used the identity 

Tr {fH^H''^) = ^ + i Z'^i^a^s) ^ P_4) 

and the anti-symmetry of the structure constant J^i'^a^s^ together with the symmetry of 

^0102 as 

The sum over the six gluon permutations of boxes is then proportional to the single color 
factor f°-f^2a3^ g^j-^^ because of Bose symmetry of the gluons, the kinematic box factor 

B''^'^'^^{quq2,q3)=B'''''\qi,q2,q3) + B'''''\q2,q3,qi)+^^^^ (D.5) 
is totally antisymmetric in the gluon indices (g^, fii), i = 1, 2, 3. 

D.l qg^ qgH and gg ggH 

The general structure of Eq. ( p.5| ) can be further restricted for the processes we are in- 
vestigating: qg^qgH and gg^ggH. In fact, in both processes, two gluons of the box 
are on-shell, and the amplitude is contracted by the two corresponding polarization vectors 
ej, while the third gluon is contracted with the conserved current 72*1 (see Eq. ( p.5|) ), in 
the qg qgH process, and with the conserved current of the two on-shell gluons in the 
three-gluon vertex, in the gg —>■ ggH case. 

This gives rise to a few simplification in the structure of Eq. (|D.5| ). In fact, a parity even, 
three-index tensor which depends on three independent momenta (here taken as the three 
outgoing gluon momenta g^, i = 1,2, 3) can be written in terms of 36 independent tensor 
structures, 9 of type g'^^^'^q^'^ and permutations, plus 27 tensors of type qi^qj^q^^- However, 
any terms proportional to gf^ g2^, or g3^ vanish by virtue of the transversity of the gluon 
polarization vectors ef* and because the current on the vertex is conserved. This leaves 
us with 14 possible tensor structures, that can be further reduced to three once we impose 
the total antisymmetry in the gluon indices {qi, fii) 

- g^'^^^qi^' Baiq2, qi, gs) - g'^'^q'^' Sa(gi, ga, gs) - g^'^'q!;' Ba{q,, g2, gi) 
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+ Qs'Qs^Qi' BbiQi^ Q2, gs) + gfgf Bb{q2, gs, ^i) + gfgf # B^igs, qi, ^2) 

- q^'q^'q^' B,{q,, q,, q,) - q^^q^'q^' B,{q,, q,, q^) - qi^^q^^qt' B,{q,, g^, qi) 

+ q'2'q'^q'i' B^iqu q2, q-,) - qt^'q'^'q^' 3^2^ qu qs). (D.6) 

Note that Bose symmetry implies that Be must be invariant under cychc permutations of 
its arguments 

Bc{qi, q2, gs) = Be{q2, ga, gi) = ^^(gs, gi, g2). (D.7) 

However, a convenient choice of gauge will remove the Be terms altogether, as shown in 
Sec. |3]^, using the polarization vectors defined in Eqs. ( |3.27|) and ( p.28|) . 



The scalar functions appearing in Eq. ( p.6| ) are given by 
1 



^a(gi,g2,g3) 



^q2 ■ ga 
- gi ■ g2 



Do{qi, g2, gs) + Do{q2, gs, gi) + Do{q3, gi, g2) 
^i3(g2, ga, gi) + ^i2(g3, gi, g2) - ^lalgs, g2, gO 



-4 



^3i3(g2, gs, gi) + Dsuiqs, gi, g2) - Dsisiqs, g2, gi) 



Bb{qi,q2,q3) = ^i3(gi,g2,g3) +^i2(g2,g3,gi) - ^i3(g2,gi,g3) 



+ 4 



^37(gi, g2, gs) + ^23(gi, g2, gs) + ^S8(g2, gs, gi) 



Bc{qi,q2,q3) 



+ ^26(g2, gs, gi) - ^39(g2, gi, gs) - A3(g2, gi, gs) 
[^o(gi, g2, gs) + ^o(g2, gs, gi) + ^o(gs, gi, g2) 
+ 4 L'26(gi, g2, gs) + ^26(g2, gs, gi) + ^26(gs, gi, g2 



+ ^3io(gi, g2, gs) + ^sio(g2, gs, gi) + ^sio(gs, gi, g2) 



Co(gi,g2 + gs) 

(D.8) 



(D.9) 



(D.IO) 



Further simplifications appear in the evaluation of qg—>-qgH matrix elements. In fact, 
since the polarization vectors for the two on-shell gluons are either proportional to or y'^ 
(see Eqs. (|3.27| ) and (|3.28|) ), we only need the contractions B^^f^^ = x^^x^^B^^^"^ ^.j^{qi, q2, gs), 
— Xf^-i^y^2B^^'^^ ^3, etc. The fi^, index will be contracted with the fermion current J21 
(see Eq. ( p.20|) ). Since qi, q2, gs and y span Minkowski space, J21 can be expanded as 

1 



21 



where 



detQg 



gi ■ J21 + q2- J21 + y- J21 y^ 



q2 -qsx^ + 



gi ■ gs a;'' + 



detQg 

g2 

gi ■ g2 

detQg ^ 



gi ■ g2 



gr 



(D.ll) 

(D.12) 
(D.13) 
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and they satisfy the orthogonahty relations 



u- qi = det Qs 
■ gi = , 



M ■ ga = , 

V ■ q2 = det Qs , 



M ■ = 
■ gs = 



M ■ ?/ = 

f ■ ?/ = 



(D.14) 



by virtue of the two on-shell conditions gi = and g| = 0. Note that there is no gs 
contraction in Eq. (|D.11| ), since ■ J21 = by current conservation. 

The orthogonahty of to aU gluon momenta g^, i = 1,2,3, imphes that ah contractions 
of 5^1/^2/^3 ^itJi an odd number of y wiU vanish. This leaves us with six non-zero contractions 
of the tensor box integrals, Byyu, Byy^, By^y, B^yy, B^xu and B^xv Via the Bose symmetry of 
the tensor integral in Eq. (p.6|) , the first four and the last two are related by a permutation 
of gluon momenta: 
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91,92) 
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Baiq3,q2,qi) 



(D.15) 
(D.16) 

(D.17) 
(D.18) 



(detQg)^ \ Ba (gi, g2, gs) - 

I qi-q2 
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(detQg) 
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Ii ■ 92 J 



^13(92, 93, 9i) + ^12(93, 9i, 92) - ^13(93, 92, 9i) 



29i ■ 92 92 ■ 93 



^0(91, 92, 93) + ^0(93, 9i, 92) - 92 ■ 93 92 ■ ( 93 + y ) ^0(92, 93, 9i) 



91 



+ 4gi - gs £'313(92,93,91) +£312(93,91,92) -£313(93,92,91) 
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-4 



92 ■ 93 93 ■ (91 - 92) - 91 9i ■ 92] [£13(91, 92, 93) + £12(92, g3, 9i) - £13(92, 9i, 93) 
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detQs r^^^(^^^^ £23(91, 92, 93) + £33(92, 93, 9i) + £25(92, 93, 9i) 



9i ■ 92 



£39(92, 9i, 93) - £23(92, 9i, 93) +91-92 Co{qi, ga + gs 



(D.19) 



-(detQ3)N5a(g2,gi,g3) + 



9i ■ 93 



9i ■ 92 



5a(93, 91,92) - Ba{q3,q2,qi] 
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[qi ■ q2Y J 
= -b^Uqi ^ Q2) . (D.20) 

Note that in the replacement B^^^ = —B^xuili ^ Q2), the Gram determinant, detQs, is to 
be treated as totally symmetric under interchange of gluon momenta qi, q2, and ^3. 

E Tensor and color structure of pentagons 
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Figure 13: Two five-point functions connected by charge conjugation. 



The two generic five-point functions connected by charge conjugation depicted in Fig. [13 
have the following expression 

d^k ^ f 1 .. 1 .. 1 
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— , (E.l) 

where 91, 92, 93 and are the outgoing momenta {qij = qi + qj and similar ones for qiji and 
qijin) and where we put an overall factor l/(4mf) in front to delete a corresponding term 
coming from the trace over the top quark. From charge conjugation we have 



pMiM2M3M4^^^^ g2, 93, 94) = P^^'^^^^'^^(gi, g2, 93, 94) 



(E.2) 



A I9i, 92, 93, 94J 

Finally, the color structure of the sum of the two diagrams depicted in Fig. |13| is 

Tr {en'^H^'H-') i'^2'^3M4(^^^ ^ Tr (^^t'^^t^^t"^) pf ^2'^3M4(^^^ 

= [Tr {en^H^H''^) + Tr (r^^^r^^'^a)] pMiM2M3M4(^^^ _ (E 3) 
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Further details about the color structure of pentagons are given in Sec. O 
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